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ABSTRACT 

The spectra of O and B supergiants are known to be affected by a significant form of extra 
line broadening (usually referred to as macroturbulence) in addition to that produced by stellar 
rotation. Recent analyses of high resolution spectra have shown that the interpretation of this 
line broadening as a consequence of large scale turbulent motions would imply highly supersonic 
velocity fields in photospheric regions, making this scenario quite improbable. Stellar oscillations 
have been proposed as a likely alternative explanation. As part of a long term observational 
project, we are investigating the macroturbulent broadening in O and B supergiants and its 
possible connection with spectroscopic variability phenomena and stellar oscillations. In this 
letter, we present the first encouraging results of our project, namely firm observational evidence 
for a strong correlation between the extra broadening and photospheric line-profile variations in 
a sample of 13 supergiants with spectral types ranging from 09.5 to B8. 
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Subject headings: stars: early-type — stars: atmospheres — stars: oscillations — stars: rotation — 
supergiants 



1. Introduction 

The presence of an important extra line broad- 
ening mechanism (in addition to the rotational 
broadening and usually called macroturbulence) 
affecting the spectra of O and B supergiants 
(Sgs) was initially suggested by the deficit of 
narrow lined obje c ts am o ng these types o f stars 
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Howarth et al.l 119971 ). The advent of high-quality 
spectra permitted confirmation that rotational 
broadening alone is insufficient to fit the line pro- 
files in many objects, and to investigate the possi- 
ble dise ntangling of both broadening contributions 
(see e .g. Rvans et al.ll2002l ISimon-Di'az fc Herrero 
2007). These studies definitely showed that the ef- 
fect of macroturbulence is dominant in the profiles 
of metal lines in early B Sgs. 

Although it was named macroturbulence at 
some point, the interpretation of this extra broad- 
ening as the effect of turbulent motion is quite 
improbable. The effect is present in photospheric 
lines and affects the whole profile, even wave- 
lengths close to the continuum. Therefore, what- 
ever is producing the extra broadening has to be 
deeply rooted in the stellar photosphere (and pos- 
sibly deeper) , in layers where we do not expect any 
significant velocity field in these stars. 

If interpreted as turbulent motion, macroturbu- 
lence would represent highly supersonic velocities 
in many cases (see Figure [J). This interpretation 
is incompatible with the previous statement. 

One physical mechanism suggested as the ori- 
gin of this extra broadening relates to oscillations. 
Many OB Sgs are known to show photometri c and 
spectr oscopic variability 



Lucy 



Based on this 
(|l976f ) postulated that this variability might be a 
pulsation phenomenon, and macroturbulence may 
be identified with the surface motions generated 
by the superposition of numerous non-rad ial os- 
cillations. More recently, lAerts et"aL I (l2009h com- 
puted time series of line profiles for evolved mas- 
sive stars broadened by rotation and thousands 
of low amplitude non-radial gravity mode oscilla- 
tions and showed that the resulting profiles could 



1 Based on observations made with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway, and Sweden, in the 
Spanish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 
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Fig. 1. — Macroturbulent velocities for B Sgs with 
spectral types B0-B3, measured in recent stud- 
ies. Two different symbols are used depending 
on the type of definition assumed for the macro- 
turbulent profile. Diamonds: isotropic Gaus- 



sian ([Dufton et al.l l2006t iFraser et all l2010h: tri 



angles: radial-tangential d efinition ([Lefever et al 



2007: iMarkova et al. 2008). Characteristic values 



for the sound speed in the photospheres of these 
types of stars are also indicated as a solid line. 



mimic the observed ones. In their paper, Aerts et 
al. (2009) considered the amplitudes of the pulsa- 
tion modes to be sufficiently low to allow a lin- 
ear superposition of mode velocities to derive the 
overall pulsational velocity eigenvector, i.e., non- 
linear mode coupling or other nonlinear effects 
were ignored. We refer the reader to Sect. 2.1 in 
Aerts et al. (2009) for a detailed description of the 
simulations, and a discussion about how the col- 
lective effect of low amplitude non-radial gravity 
modes can produce the inferred supersonic level of 
macroturbulent broadening when the oscillations 
are ignored in the interpretation of spectral lines. 
Stellar oscillations are a plausible explanation for 
the extra broadening in O and B Sgs, but so far 
there is no direct observational evidence confirm- 
ing this hypothesis. 

Two years ago we began an observational 
project aimed at investigating the macroturbulent 
broadening in O and B Sgs, and its possible con- 
nection with spectroscopic variability phenomena 
and stellar oscillations. In this letter, we present 
the first encouraging results, showing observa- 
tional evidence for a strong correlation between 
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this extra broadening and photospheric line pro- 
file variations in a sample of 13 Sgs with spectral 
types ranging from 09 to B8. 

The letter is structured as follows: the sample 
of stars and the spectroscopic observations used 
for this study are presented in Sect. [2J we char- 
acterize the line-broadening and the line-profile 
variations of photospheric lines in Sects. [3] and 
21 respectively; a connection between the extra 
broadening and line profile variations is presented 
in Sect. [5j Finally, we discuss the possible relation 
to stellar pulsations in Sect. [6] 

2. Observational data set 

We obtained time series spectra of a selected 
sample of 11 bright late-0 and early-B Sgs during 
two observing runs (2008/11/05-08, 2009/11/09- 
12) with the FIES cross-dispersed high resolution 
echelle spectrograph on the 2.5m Nordic Optical 
Telescope at Roque de los Muchachos Observatory 
on La Palma (Canary Islands, Spain). We used 
FIES in the medium resolution mode (R=46000, 
(5A=0.03 A/pix). The sample was complemented 
with two OB dwarfs and two late B Sgs. The ex- 
posure times were chosen so as to reach at least 
SNR=200 (measured in the range 4500-4600 A). 
The list of observed stars, along with their spec- 
tral classification, number of spectra obtained for 
each target, time span of the corresponding time- 
series, and observing run in which they were ob- 
served are indicated in the first columns of Table 
[TJ Note that four of these stars were observed in 
both campaigns. 

The spectra were reduced with the FIEStooH 
package in advanced mode. A proper set of 
bias, flat and arc frames obtained each night 
was used to this aim. The FIEStool pipeline 
provided wavelength-calibrated, blaze-corrected, 
order-merged spectra of high quality. Next, the 
barycentric correction was performed. We used 
the interstellar Na I D doublet at 5890 A to check 
the accuracy in this correction, and found an 
agreement in the line doublet position for all the 
time series spectra of each star better than 0.5 
kms -1 . 

Once the lines of interest for this study were 
identified, a local normalization of selected regions 



of the spectra was performed. For each line, the lo- 
cal continuum used for the normalization was cal- 
culated from a linear fit between points inside two 
continuum regions adjacent to the line. The same 
continuum windows (selected manually in the av- 
erage spectrum) were used for all the spectra of a 
given star. 

3. Characterization of the line broadening 
in photospheric lines 

Several types of broadening mechanisms are 
traditionally considered to produce the final ob- 
served line profiles of photospheric metal lines 
in massive stars: instrumental, natural, thermal, 
microturbulent, macroturbulent and rotational 
broadening. In the case of high resolution^ spectra 
of late-0 and early B Sgs, rotational and macro- 
turbulent constributions dominate. 



W e used the Fourier transform technique (|Grav 
1976: see also ISimon-Diaz fc Herrerol 12007 1 for a 
recent application to the spectra of O and B stars) 
to disentangle the contributions from rotational 
and macroturbulent broadening. The results of 
the analysis are presented in Table [TJ We per- 
formed the analysis for each of the time series spec- 
tra, obtaining v sin i as indicated by the first zero 
of the Fourier transform. The range and median 
of derived v sin i values are indicated in columns 6 
and 7 of Table Q] Note that the dispersion in the 
obtained usinz values is between 10% and 30%, 
depending on the star. Whether this dispersion 
is real or an effect of nois e is not clear from this 
data s et. As outlined by ISimon-Diaz fc Herrero 
(|2007l ). the correct identification of the first zero is 
complicated in cases of low SNR and for a large ex- 
tra broadening contribution. This will be explored 
in more detail in a future investigation. 

Next, we applie d a goodness-of-fit method (viz. 
Ryans et al.ll2002l ) to quantify the contribution of 
the extra broadening in each spectrum from the 
time series. We considered as free parameters the 
equivalent width, the radial velocity, wsini, and 
the parameter defining the extra broadening. We 
allowed v sin i to vary in the fitting process in or- 
der to compare with results from the FT anal- 



2 http: / / www.not.iac.es / instruments / fies / fiestool/FIEStool.html 



3 The instrumental broadening contribution can be mini- 
mized by using high resolution spectra. In our FIESONOT 
data set (R=46000) the instrumental Gaussian profile has 
a FWHM ~6.5 kms' 1 . 
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Table 1: Columns 1-5: Targets, spectral classification, observing campaign, time span of the observations, 
and number of spectra. Columns 6-11: Line-broadening and LPVs characterization in the Si III 4567 or 
O III 5592 lines for our stellar sample. LPVs characterized with the peak-to-peak amplitudes of the first and 
third moments of the lines. AT in days; vsini, Ort, { v )ppi and (^ 3 }pp 3 hi kms -1 . Uncertainties in (v) pp , 
and (u 3 )p P 3 are ^0.10-0.30 and ^0.4-1.2 kms -1 , respectively 



■usina(FT) x 2 fitting LPVs 



Star 
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HD 207198 


09 Ib-II 


FI09 


3.14 


14 


59- 


-64 


62±2 


50±2 


100±0 


3.41 


26.9 


HD 209975 


09.5Iab 


FI08 


3.04 


21 


56 


02 


59±2 


50±6 


95±5 


10.13 


46.3 






FI09 


3.18 


18 


52- 


-58 


57±2 


51±4 


100±5 


11.65 


45.4 


HD 37742 


09.71b 


FI09 


3.19 


30 


90 


140 


117±10 


115±15 


117±13 


12.01 


62.8 


HD 204172 


B0 lb 


FI09 


3.14 


14 


38- 


-68 


58±11 


58±6 


85±5 


11.80 


44.0 


HD 37128 


B0 la 


FI08 


3.01 


60 


45- 


-64 


55±9 


57±8 


90±6 


12.30 


43.0 






FI09 


3.18 


28 


32 


60 


49±8 


54±9 


90±9 


13.02 


46.5 


HD 38771 


BO. 5 la 


FI08 


2.99 


48 


42 


-57 


50±6 


48±5 


90±6 


9.38 


36.6 






FI09 


3.17 


25 


30 


-55 


48±8 


52±7 


90±5 


10.56 


39.1 


HD 2905 


BC0.7Ia 


FI08 


3.07 


30 


44 
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90±3 


10.07 


39.8 
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39.6 


HD 24398 


Bl lab: 


FI09 


3.19 


34 


27 


43 


34±5 


38±2 


65±0 


3.15 


25.7 


HD 190603 


B1.5Ia+ 


FI08 


1.96 


17 


30 


47 


46±5 


49±4 


55±3 


2.18 


24.6 


HD 14818 


B2Ia 


FI08 


3.00 


16 


32 


47 


41±6 


41±5 


70±3 


7.84 


37.1 


HD 206165 


B2Ib 


FI09 


3.14 


16 


29 


41 


38±2 


38±2 


60±0 


5.63 


28.2 


HD 191243 


B5Ib 


FI09 


2.10 


9 


15 


25 


18±3 


18±2 


30±0 


1.36 


13.4 


HD 34085 


B8Iab: 


FI09 


3.08 


21 


14 


33 


23±4 


29±4 


35±5 


6.18 


19.8 


HD 214680 


09 V 


FI08 


3.12 


31 


13 


23 


18±3 


16±2 


40±5 


2.12 


17.4 


HD 37042 


B0.5V 


FI08 


2.92 


20 


32 


34 


33±1 


33±0 


15±5 


0.84 


10.4 



ysis. Two possibilities for the characterization of 
the extra broadening were considered: an isotropic 
Gaussian profil e, and a Ga ussian radial-tangential 
one (6rt, see lGravlll976l ). In the latter, we as- 
sumed that the radial and tangential components 
provide equal contributions to the final profile. 
The analysis using an isotropic Gaussian profile 
for the macroturbulent broadening resulted in (a) 
a large dispersion in the vsini values obtained 
from the various spectra in each time series; (b) 
systematically lower v sin i values than those de- 
termined through the FT method (a factor of 0.5 
in many cases, sometimes even "zero" rotation) 
for all stars with dominant extra broadening. On 
the other hand, when assuming a Gaussian radial- 
tangential profile, we found that (a) the disper- 
sion in the vsini and Ort values was rather low; 
(b) differences between the v sin i values obtained 
from the goodness-of-fit method and using the FT 
technique are < ±20%. We thus conclude that a 
Gaussian radial-tangential profile is better suited 
to characterizing the extra line broadening than 
an isotropic Gaussian profileOl The corresponding 



mean and standard deviation for v sin i and Ort as 
obtained from the x 2 fit are indicated in columns 
8 and 9 in Table Q] Note the good agreement in 
the derived values for the four stars observed in 
both 2008 and 2009. 



Sim ilarly to previous studies (e.g 



2006: Lefever et al. 2007: Maxkova et all 



Dufton et al 
2008; 



Fraser et al.ll2010l ). we found that the size of the 
extra broadening is generally larger than the ro- 
tational contribution in all the Sgs. This is not 
the case for the B0.5V star HD 37042, where the 
total broadening is mainly produced by the ef- 
fect of the stellar rotation (note, however, that 
a certain extra broadening is also needed in this 
case). The other dwarf star, HD 214680, is a spe- 
cial case, since it has very low vsini. For such 
a low vsini values, microturbulence produces a 
significant contribution to the total broadening, 
which then is included in the measured extra 
broadening. Note that this argument may also 
explain the extra broadening needed for HD 37042 
(and the other two late-B Sgs with low vsini). 



subtle. 



4 We also found that the former generally produces a better 
final fit to the observed profiles, though the differences are 
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Fig. 2. — Representative examples for the variation of the first moment (radial velocity placed at average 
zero) and third moment (skewness of the line) for some of the stars in the studied sample derived from the 
Si III 4567 or O III 5592 line profiles. Horizontal lines in upper plots show the accuracy in radial velocity 
associated with the instrumental setting used for FIES@NOT. 



4. Characterization of line profile varia- 
tions in photospheric lines 

In agreement with earlier studies of spectro- 
scopic variability in O and B Sgs (e.s. lEbbets 



1996 : IPrinia et al.ll 
2004 ; iKaufer et al 



996, 



19821: iHowarth et al.ll993tlFullerton. Gies. fc B olton 



20041 120061 fMorel et al 



2006t iMarkova et al 



20051 



2008), we found clear signatures of line profile 
variations (LPVs) for all the Sgs considered in 
our sample. To investigate these LPVs quan- 
titatively we computed the first, (u), and third, 
(f 3 ), nor malized velocity mo ments (for definitions, 



see, e.g.. lAerts et al.ll2010aD from t he Si III 456 7 
or O III 5592 lines with FAMIAS^I (jZimal 120081 ). 
These moments are related to the centroid veloc- 
ity and the skewness of the line profile, respec- 
tively, and are well suited for deciding whether an 
observed line profile is subject to time-dependent 
line asymmetry (note that their values are zero for 
purely symmetric profiles). 

Four representative examples of the tempo- 
ral behavior of the first and third velocity mo- 
ments are presented in Fig.[5J The associated un- 
certainties (not included in the plot) are — 0.1- 
0.4 kms -1 , and —1000-4000 km 3 s" 3 , respec- 
tively. In the case of the dwarf star HD 37042, 



5 A software package developed in the framework of 
the FP6 European Coordination Action HELAS 



<http://' 



www.helas-eu 



■org/). 



the (v) values are fairly constant and close to 
zero. The maximum dispersion in velocity for 
this star is —1 kms , of the order of the ac- 
curacy associated with the instrumental setting 
used for the FIES@NOT observations (indicated 
as dashed horizontal lines). All the other stars 
in Figure [2] show (v) variations above this sig- 
nificance level, with maximum peak-to-peak am- 
plitudes reaching ~1Q— 12 kms -1 in some of the 
cases. The minimum variations are found for 
HD 214680 (the other luminosity class V object 
considered in this study), HD 191243 (B5Ib) and 
HD 190603 (B1.5Ia + ), with maximum amplitudes 
slightly larger than the significance level. Curi- 
ously, the star HD 190603, expected to be an ex- 
treme object from its spectral classification, is one 
of the cases with smaller variations. We indicate 
in Table Q] (columns 10 and 11) the peak-to-peak 
amplitude of the first and third moment variations 
measured for each of the considered targets. 

5. The "macroturbulence"-LPV connec- 
tion 

We investigate the possible connection between 
macroturbulent broadening and LPVs in our sam- 
ple of stars. Such a connection should appear 
in case macroturbulent broadening is produced 
by any time-dependent ph ysical phenomena (e.g. 
non-radial oscillations, see lAerts et al. 2009). In 
Fig. El (upper pannels), we plot the average size of 
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the macroturbulent broadening, (0rt), versus the 
peak-to-peak amplitude of (v) and (v 3 ) variations 
for each of the stars studied. Results from the four 
stars observed in both campaigns are conected 
with solid lines. A clear positive correlation is 
present in both cases. To our knowledge, this is 
the first clear observational evidence for a connec- 
tion between extra broadening and LPVs in B and 
late O Sgs. Particularly remarkable is the Ort - 
(v 3 ) correlation: the larger the extra broadening, 
the more asymmetric line profiles can be found in 
the time series. Note that this does not mean that 
lines with a significant macroturbulence contribu- 
tion are always asymmetric since (i> 3 ) oscillates 
between positive and negative values over time. 

We also present in Fig. [3] (bottom p annels) sim- 
ilar plots with data from Table 1 in Aerts et al 



(|2009h , based on simulations of line profiles broad- 
ened by rotation and by hundreds of low amplitude 
non-radial gravity mode pulsations. These simu- 
lations consider various combinations of the incli- 
nation angle («), the projected rotational velocity 
(wsini), and the i ntrinsic amplitude of the modes 
denoted by a (see Aerts et all 120091 for a defini- 
tion). We present the maximum values of v mac 
obtained for each set of simulated profiles. Similar 
trends for i> mac as in Fig. [3] are found using aver- 
age or minimum values for the macroturbulence 
parameter, but with a different (smaller) scale in 
the y-axis. 

The simulations lead to clear trends which are 
compatible with spectroscopic observations. 

6. Is macroturbulent broadening in OB- 
Sgs caused by pulsations? 

Non-radial oscillations have been often sug- 
gested as the origin of LPVs and photosperic lines 
in OB Sgs, as well as the driver of large scale wind 
structures (see references in Sect.HJ; however, a 
firm confirmation (by means of a rigorous seis- 
mic analysis) has not been achieved yet. From a 
theoretical point of view, g-mode oscillations were 
not initially expected in B Sgs because the radia- 
tive dampi ng in the core wa s suspected as being 
too strong. ISaio et aL ( 20061 ) detected simultane- 
ous p- and g-modes in HD 163899 (B2Ib/II) using 
data from the MOST satellite. These authors also 
computed new models showing that g-modes can 
be excited in massive post-main sequence stars, as 




024681012140 20 40 60 

<v> pp (km s -1 ) < v3 >pp' /3 ( km s "') 

Fig. 3. — (Top) Empirical relations between 
the average size of the macroturbulent broad- 
ening ((Ort)) a n d the peak-to-peak amplitude 
of the first and third moments of the line- 
profile. (Bottom) S imilar relations are found from 
the simulations bv lAerts et al. (|2009h . 



the g-modes are reflected at the co nvective zone 



associ ated with the H-burning shell. iLefever et al 



(|2007l ) presented observational evidence of g-modc 
instabilities in a sample of photometrically vari- 
able B Sgs from the location of the stars in the 
(log 5, logT cff )-diagram. 

These results, along with our observational con- 
firmation of a tight connection between macrotur- 
bulence and parameters describing observed LPVs 
render stellar oscillations the most probable phys- 
ical origin of macroturbulent broadening in B Sgs; 
however, it is too premature to consider them as 
the only physical phenomenon to explain the un- 
known broadening. 

Recent studies of high resolution spectra of OB 
stars indicate that the appearance of a signifi- 
cant macroturbulent broadening contribution is 
not only limited to B Sgs but occurs also in O 



stars of all luminosity clas ses (|Sim6n-Diaz et al 



20ld; iMarkova et~al] l2010h as well as 



m mam- 



sequence B-type pulsators (e.g. Morel et "al1 l2006') . 
Moreover, with the recent discovery of a strange- 
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mode oscillation in the B6Ia supergiant HD 46005 
( Aerts et al. 2010bl ) and of stochastic oscillations 
in the 08V star HD 46149 (iDeeroote et al.ll2010h 
and t he j5 Cephei star V1449Aql (jBelkacem et all 
2009) from high-precision space photometry, ob- 
servational evidence of atmosphere phenomena 
due to oscillations, which must cause some level 
of pulsational broadening, is growing. 

Future investigations of the correlation between 
LPVs and macroturbulcnce, and its connection to 
stellar oscillations in the whole realm of massive 
stars, will be key to establishing the role played 
by the different physical processes contributing to 
the spectral line formation in these types of stars. 
Different models including various combinations of 
pulsational, rotational and wind effects should be 
considered, keeping in mind that an appropriate 
physical mechanism must be able to explain the 
observed correlation presented in this study. 
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